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ABSTRACT 

Several relatively bright, persistent X-ray sources display regular pulses, with periods in the 
range of (0.7 - 10) x 10^ s. These sources are identified with massive close binaries in which a 
neutron star accretes material onto its surface. The observed pulsations in all of them, but one, 
are unambiguously associated with the spin period of the neutron star. Analyzing possible 
history of these pulsars I conclude that the neutron stars in these systems undergo spherical 
accretion and their evolutionary tracks in a previous epoch contained three instead of two 
states, namely, ejector, supersonic propeller, and subsonic propeller. An assumption about a 
supercritical value of the initial magnetic field of the neutron stars within this scenario is not 
necessary. Furthermore, I show that the scenario in which the neutron star in 2S 01 14+650 is 
assumed to be a magnetar descendant encounters major difficulties in explaining the evolution 
of the massive companion. An alternative interpretation of the spin evolution of the neutron 
star in this system is presented and the problem raised by association of the lO'* s pulsations 
with the spin period of the neutron star is briefly discussed. 
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1 INTRODUCTION 

A newly formed neutron star is presumed to rotate rapidly with 
a period of a fraction of a second. Its rotational rate then de- 
creases, initially by the conventional spin-powered pulsar energy- 
loss mechanism (ejector state), and later by means of the interac- 
tion between the magnetosphere of the neutron star and the stellar 
wind of its companion (propeller state). As the spin period of the 
neutron star reaches a critical value the accretion of material onto 
its surface starts (accretor state) a nd the star switches on as an X- 
ray pulsar (for a review see, e.g., Bhattacharva & van den Heuvell 



ll99lHlben. Tutukov & YungelsoiJ 



nil). 



This scenario has been numerically studied by Urpin, Ko- 
nenkov & Geppert (1998). Their calculations have been performed 
under assumptions that the initial magnetic field of the neutron star 
is subcritical (i.e. Bo < B„ = mlc^/eTi = 4.4 x 10''' G) and has 
a crustal origin. They have also assumed that the accretion onto 
the stellar surface starts as soon as the centrifugal barrier at its 
magnetospheric boundary becomes ineffective. The period at which 
this condition is satisfied is defined by equating the corotational ra- 

dius of the star, rc = [GMP'^IAn \ , to its magnetospheric radius 
=K{frim i2GMf'\ that gives 
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(for a discussion see, e.g., Pringle & ReesI 1 19721 : 
llllarionov & Sunvaevlll975l : IStella. White & Rosned ligseV Here Ze- 
is the parameter accounting f or the geometry of the accretion flow, 
which ranges from 0.5 to 1 jOhosh & Lamblll978l) . and tsd is the 
total spin-down time of the star in the ejector and propeller states. 
M, Ps, and fi are the mass, spin period, and the dipole magnetic 
moment of the neutron star, respectively. TO = nr^p„V^^i is the 
mass with which a neutron star interacts in a unit time as it moves 
through the wind of a density p„ with a velocity VkI = V^ns + ^w' 
where V„s is the linear velocity of the neutron star's orbital motion 
and V„ is the wind velocity of its massive companion. The param- 
eter ro = 2GM/V^^| denotes the Bondi radius of the neutron star, 
and m = M/IAMq. The acceleration torque applied to the star 
in the accretor state has been evaluated as K^u = l^JiiGMr^y^'^, 
where ( is the efficiency parameter, which in the calculations has 
been normalized to 0. 1 . Finally, the evolutionary tracks of neutron 
stars have been computed on a time scale ~ 2 x 10' yr, which 
corresponds to the life-time of a SOMp star on the main sequence 
tehattacharva & van den Heuvellll991l ). 

According to their results the spin period of a neutron star 
rapidly increases during the ejector and propeller states to a max- 
imum value ~ Pcd(tsd) and gradually decreases as the accretion of 
material onto their surface starts (see Figs. 2 and 4 in lUrpin et al.l 
1998). The maximum period, which a neutron star with Bq < 
is able to reach on a time scale tsd < fms under the above condi- 
tions is limited to a few hundred seconds. In particular, for the case 
Bo = 10" G, the longest period is ~ 500 s, which the neutron star 
reaches at TO ~ 10"'^ Mq yr"'. The magnetic field of the neutron 
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Table 1. Very long-period X-ray pulsars 



Name 


Sp. type 




Po,b, d 


log Lx* 


Ref." 


J1700()6-4157 




715 




34.7 


[1] 


0352+309 


BO Ve 


837 


250 


34.7 - 35.5 


[2, 3] 


(X Per) 


(09.5 Ille) 










J1037.5-5647+ 


BO V-IIIe 


862 




34-35 


[4,5] 


J2239.3+6116 


BO Ve 


1247 


262.2 


33-36 


[6, 7, 8] 




(B2 nie) 










J0103.6-7201 


05 Ve 


1323 




35.3-36.8 


[9] 


J0146.9+6121 


Bl Ill-Ve 


1412 




34.6-36 


[10] 


(V831Cas) 












2S 01 14+650 


Bl la 


10008* 


11.6 


35.7-36 


[11, 12] 



* Lx is expressed in erg s ' , 

" Ref erences to the tabl e: [1] iTori i et al.' Iwf); [2] 'Pelgado-Marti et alJ 
j200l[) : [3] Haberl et al ]h998l) : [4] Reig & Roche ( 1999); [5] Motch e Tal] 
'l997); [6] in 't Zand et al. (2000); [7] in 't Zand et al. (2001^ 
[8] Ziolkowski (2002); [9] Haberl & Pietsch ( 2005.); [101 .Haberl & Sasaki, 
i2OO0i) ; llll [ReigetalJjl996l) ; 112l lHall et al] j2000l) . 

* The Be/X-ray transient pulsar. 

* An association of the 10008 s pulsations with the spin period of th e 
neutron star has recently been challenged by iKoenigsberger et aP j2006l) . 
see text. 

star at the end of the propeller state is a factor of 3 weaker that its 
initial value, B(t^) ~ 0.3So, and decreases by almost an order of 
magnitude on a time scale of binary evolution, B{t„^s) ~ 0. ISo- 

The long-period pulsars within this scenario are expected to 
be relatively weak X-ray sources. Indeed, solving Eq. l[T) for SR 
and putting the result to Lx = djiGM/r„s one finds 

\-7/3 , 



ix <4x 



(sooi) (a3^) (to^) 



(2) 



where r„s is the neutron star radius. The pulsars remain in the 
low-luminosity state as long as the mass transfer is operated via 
the wind-fed mechanism. As the massive companion overflows its 
Roche lobe the mass-transfer rate increases. However, this transi- 
tion will unavoidably be accompanied with a formation of an ac- 
cretion disc and, the refore, a rapid spin-up of the neutron star o n a 
time scale of (Davidson & Ostrikedll973l ; ]Ghosh & Lambll 19781) 



Af- 



IttI 



(3) 



where / is the neutron star moment of inertia. Under the conditions 
of interest, i.e. / = lO"*^ Ins g cm^ and TO =10'^ SjJis g s"' , the spin- 
up time scale is < 1000 yr. This indicates that the probability to 
detect the pulsar in this stage is almost negligible. 

There are, however, seven X-ray pulsars whose periods sub- 
stantially exceed 500 s (see Tab.[T]l. These sources are associated 
with massive close binaries in which a magnetized neutron star ac- 
cretes material onto its surface. All of them but one are persistent 
X-ray emitters with an average luminosity > 10'* ergs"', which 
implies 211 > lO^'-MQyr"'. Observations of a cyclotron line in 
the X-ray spectrum of X P er suggest that the surface field of the 
neutron star is 3 .3 x 10'^ G dCobum et alJlOOll) . There is some evi- 
dence (see, e.g.. lLi & van den Heuvell 19991) that the magnetic field 
of neutron stars in other sources listed in Tab.[T]is of the same order 
of magnitude. 

The origin of these sources within the above scenario is rather 
unclear. Indeed, the value of P^dihii) for 6(4d) - 0.3ficr and 
5R > 10"'- Mq yr"' is limited to < 300 s. For a neutron star with 
Bo < to reach a period of ~ lO' s the strength of the wind should 



be < 6 X lO"'"* Mq yr"' . In this case, however, the luminosity of the 
pulsar during the wind-fed accretor state is < 7 x lO'^ergs"', i.e. 
by a factor of 130 smaller than the average X-ray luminosity of 
the sources listed in Tab.[T] On the other hand, a possibility that 
the massive star overflows its Roche lobe at least in the case of 
X Per, J2239.3+6116, and 2S 01 14+650 can be excluded. Hence, 
the only way to explain the origin of the long-period pulsars within 
the above scenario is to invoke a rather controversial assumption 
that the mass-loss rate of their massive companions has recently 
been increased by more than two orders of magnitude. Finally, the 
lO'' s spin period of neutron star in 2S 01 14+650 represents an ex- 
ceptional case which within the above scenario appears to be a mir- 
acle. 

A possibility to improve the situation by assuming that the 
neutron star i n 2S 0114+650 was bom as a magnetar (i.e. So » 
ficr. see, e. g., 'Thompson & Duncanlll993l) has been discussed by 
iLi & van"d en Heuvel ( 1999). They have shown that the spin-down 
time scale of a star with Bq ~ 10" G to a period ~ lO'' s under the 
condition TO <, 10'** gs"' is close to 10' yr, which is comparable 
with the character istic time of magnetic field decay of magnetars 
dColpi etal] 12000 ^. Following this finding they have suggested a 
scenario in which the neutron star in 2S 01 14+650 is assumed to 
be a magnetar descendant undergoing a spherical accretion with a 
very low angular momentum transfer rate. They have also pointed 
out that this scenario, being applied to other long-period pulsars, 
may shed a new light to the origin of strong magnetic field of 
neutron stars in such systems as A0535+26, VelaX-1, GXl+4, 
4U 1907+09, 4U 1538-52, and GX 301-2. 

It appears, however, that the rotational rate of a neutron star, 
which forms during the first supernova explosion in a massive 
binary system, is insufficient f or its magnetic field to be ampli- 
fied over the critical value ffhompson & Murravl200ll ; iHeger et"ail 
l2003l : |Petrovic et al.ll200^ . It is more likely, that magnetars form 
at the latest stages of binary evolution, namely, during the second 
supernova explosion or, the most probably, during a coalescence 
of two neutron stars (see, e.g., [Price & Rosswog|l2006 „ and refer- 
ences therein). In this case, however, the fraction of magnetars in 
binary systems can unlikely exceed 1 per cent and most of them 
are expected to have a black hole companion (for a discussion see, 
e.g., Popov & Prokhorov 2006). Moreover, the kicks which mag- 
netars are expected to get dur ing their formation a re too large for 
the binary system to survive (" Wheeler et i^bOOOl. and references 
therein). In this light, a probability for a magnetar to be accompa- 
nied with a massive main sequence star appears to be almost negli- 
gible. 

While these are not compelling arguments against the long- 
period pulsars to be descendants of magnetars, they suggest that 
alternative possibilities to solve the problem might be more fruit- 
ful. One of them is discussed in this paper. Namely, I show 
that the spin periods of long-period pulsars can be expla i ned in 
terms of evolutionary tracks constructed by iDavies etd] h979l) 

I n [ 1 I 1 

and Davies & Pringle ( 1981) and recently improved by Ikhsanov| 
(2001a) provided the neutron stars in these systems undergo spher- 
ical accretion. These tracks contain an additional, the so called sub- 
sonic propeller, state, which has not bee n tak en into account in the 
models considered bv lUrpin etal] ( ll998h and lLi & van den Heuvell 
( Il999l) . The spin period of a neutron star during this state increases 
in a relatively short time, Tj <c 4^, to a value expressed by Eq. j24b . 
which under the conditions of interest lies in the interval 10' - lO* s. 
The assumption about the supercritical value of the initial magnetic 
field of neutron stars for the interpretation of long-period pulsars 
within this scenario is not required. 
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2 ACCRETION FLOW GEOMETRY 



The state of neutron stars in the long-period pulsars is unambigu- 
ously identified with an accretor. The equation governing the spin 
evolution of stars in this state reads jPavidson & Ostrikej [19731 : 
iGhosh&Lambll 19781) 



d 1 
at 



(4) 



where K^^ and A'sd are the acceleration and deceleration torques ap- 
plied to the neutron star, respectively. The deceleration torque is 
associated with interaction between the stellar magnetic field and 
material situated at a distance r ^ r^. The average value of this 
torque does not significantly depend on the geometry of the accre- 
tion flow bey ond the magnetosphere and in t he general case can be 
expre ssed as jLvnden-Bell & Pringlelll974 IWanlll98lh llTpunovl 
Il992h 



(5) 



where (kt < 1) is the dimensionless parameter of the order of unity. 
In contrast, the acceleration torque is very sensitive to the flow ge- 
ometry and varies from its maximum value, A'^„ ^ SR s/GMr^, in 
the case of a disc, to a significantly smaller value, 

K:^^ ^^9jm„br|, (6) 



in the case of a spherical accretion toavidson & Ostrikeil [19731 : 
IWang|ll98ll) . Here n^ih = '^^IPoA, is the orbital angular veloc- 
ity and f orb is the system orbital period. The parameter ^ is the 
factor by which the angular momentum accretion rate is reduced 
due to inhomogeneities (the velocity and densit y gradients) in the 
accretion flow. Numerical simulations (see, e.g. jAnzer et '^ 1 19871 : 
iTaam & Frvxe ll 1988; Ruffert 1999) suggest that the average value 
of this parameter is < ^ >= 0.2. 

As follows from Eq. Q, the spin period of an accreting neu- 
tron star evolves to the so called equilibrium period, which is de- 
fined by setting K^a = K^d. The values of the equilibrium period in 
the case of a disc and spherical accretion, respectively, are 



P^^ ^ 18 s X /f" 



and 



psph 
^ eq 



910 s X yt,"" e 



1/2 £-1/2 „,-3/2 qif-l/: 



(3.3 X 1012g) ^ 



(8) 



UoOkms-i/ \250d/ ' 



UOOkms-i/ 

where ^250 = ^"0^/250 days, and ^0.2 = f/0.2. The normalization 
of parameters in these equations is appropriate for X Per, which is 
the best studied persistent long-period pulsar. 

As easy to see, the value of P^^ under the conditions of interest 
is much smaller than the spin period of neutron stars in the pulsars 
listed in Tab.[T] Hence, if a persistent accretion disc in these systems 
existed the acceleration torque applied to the neutron star would 
significantly exceed the deceleration torque. The neutron star in this 
case would regularly spinning-up at a rate P ^ Pl^Ji s/GMr^/lnl ~ 
-5x 10"* s s"i, which implies the lifetime of the long-period pulsars 
to be < 1000 yr (see Eq. O. In particular, if the neutron star in 
2S 01 14+650 accreted material from a disc its spin period would 
substantially increase on a time scale of only 100 yr, which is 4-5 
orders of magnit ude smaller than t he average lifetime of accretors 
in massive stars dUrpin etalJI 19981) . Thus, a probability to observe 
neutron stars accreting from a disc at a rate ~ 10'^ g s"' during their 
long-period stage is almost negligible. 



Furthermore, observations of X-ray source in X Per show 
no evidence for a regular spin-up of this pulsar. Instead, it ex- 
hibits apparent erratic pulse frequency variations on a time scale 
of a few days, which are superposed with a 10-20 y r spin-up/spin- 
down trends around the average period of 837 s ( lHaberl|[T994l: 
iDelgado-Marti et al]|200lL and references therein). The transitions 
between the spin-up and spin-down stages occur without any sig- 
nificant variations of the system X-ray luminosity, and therefore, 
cannot be associated with the transitions of the neutron star be- 
tween the accretor and propeller states. The spin-down events in 
this case could occur only if the value of acceleration torque ap- 
plied to this star is much smaller than A'^^,, which argues against the 
presence of a persistent accretion disc in this system. 

In contrast, the assumption about spherical geometry of the 
accretion flow allows us to interpret the spin period of the neutron 
star in X Per in terms of the equilibrium period, P''J^, provided the 
average (on a time scale of > 20 yr) value of the wind velocity is 
350 - 400kms"i. The spin-up/spin-down behaviour of the source 
within this scenario can be associated with apparent variations of 
the stellar wind velocity (e.g., due to activity of the massive com- 
ponent and the orbital motion of the neutron star whose trajectory 
is inclined to the plane of decretion disc of the Be-companion) and, 
possibly, the flip-flop instability of the accretion flow (see, e.g., 
iTaam & Frvxellll 19881) . The spherical accretion model is also ef- 
fective for the interpretation of other long-period pulsars provided 
their orbital periods are of the same order of magnitude as that of 
X Per, and the average wind velocity is ~ 400 - 800 km s"i . 

The only exception is 2S 0114+650. The condition P^ = P''J^ 
in the case of this source implies V^ei ^ 3000 km s"', which by a 
factor of 2-3 exceeds the upper limits to the stellar wind velocity of 
B-type sta rs inferred from the UV obs e rvations (800- 1500kms~i, 
see, e.g., iBemacca & Bianchil Il98ll : ISnowl Il98ll : iMarlborougl] 
Il982h . On the other hand, assuming that the wind velocity of the 
massive companion is limited to < 1500kms"i one finds /'^Jj'' < 
26 minutes. This indicates that either the neutron star is a very 
young accretor and its period significantly exceeds the equilib- 
rium one, or the value of the parameter f in the particular case of 
this system is significantly smaller than its average value inferred 
from the numerical simulations (see above). However, it cannot be 
excluded that this apparent discrepancy occur because of a mis- 
taken association of the lO"* s pulsati ons with the spin period o f 
the neutron star. As recently shown bv lKoenigsberger et al.l ( I2OO6I) . 
these pulsations may reflect a modulation of the B-supergiant wind 
caused by tidal interaction between the non-synchronously rotat- 
ing binary components. As the material of the wind is captured 
by the neutron star and is accreted onto its surface the X-ray lu- 
minosity of the source would be modulated with the same pe- 
riod, which under the conditions of interest lies within the inter- 
val of 2-3 hours. This scenario weakens the association of lO"* s 
pulsations with the spin period of the neutron star and, there- 
fore, leaves the questi on about the rotational rate of the star open 
(for a discussion see, iKoenigsberger et al] |2006 |). As mentioned 
above, the condition P^ = P^tq for reasonable values of Vrei pre- 
dicts Ps ~ 10 - 26 minutes. In this light, it is interesting to note 
that a detection of ~ 15 m inut es pulsations has been reported by 
IKoenigsberger et al] ( Il983h and lYamauchi et al] ( ll99Cll) . These pul- 
sations, however, have not been found in later X-ray observations 
of this system. 

Summarizing this section I can conclude that a presence of ac- 
cretion disc in the long-period pulsars is very unlikely. The neutron 
stars in these systems are, therefore, undergoing spherical accre- 
tion. 
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3 SUBSONIC PROPELLER 

Let us now address the main question, namely, what is the maxi- 
mum period to which a neutron star undergoing spherical accretion 
can be spun-down on a time scale t < /n,s ? A c omprehensive analy- 
sis of this question ha s been first presented bv lDavies et alj i 19791) 
and lOavies & Pringld ( Il98lh . As they have shown the spherically 
accreting neutron star in the state of propeller is surrounded by a 
hot quasi-stationary turbulent atmosphere. The atmosphere forms 
as soon as the pressure of relativistic wind ejected by the neu- 
tron star in the spin-powered pulsar state can no longer balance 
the pressure of the surrounding material, and the latter, penetrating 
to within the accretion radius of the star, interacts with the stellar 
magnetic field. This interaction leads to formation of the magne- 
tosphere with the equatorial radius which is defined by equat- 
ing the ram pressure of the flow with the magnetic pressure due to 
dipole field of the neutron star. The kinetic energy of the flow at 
the magnetospheric boundary is converted into its thermal energy 
in the adiabatic shock. The temperature of material in the shock in- 
creases to the adiabatic (or the so called free-fall) temperature, = 
GMiUp/k^r, where lUp and are the proton mass and Boltzmann 
constant, respectively. The heated gas expands with the velocity 
Vff = (2GM/r)''^, which significantly exceeds the sound speed in 
the material captured by the star. The gas expansion is, therefore, 
leading to a formation of the back-flowing shock, which propagates 
through the flow and heats it up to the adiabatic temperature. Un- 
der the conditio n TO < 9jJ„ ^ 2 x 10"* gs"' m (V,ei/10'* cms"') 
( Ilkhsanovll2002 ') the atmosphere is extended up to the Bondi radius 
of the star, rQ, at which the thermal pressure of the heated material 
is equal to the ram pressure of the surrounding gas, ~ PooV^^p The 
formation of the atmosphere prevents the stellar wind from pene- 
trating to within the Bondi radius of the neutron star. As the neutron 
star moves through the wind of its companion the surrounding gas 
overflow the outer edge of the atmosphere with the rate TO and the 
mass of the atmosphere is conserved. 

As long as the magnetospheric radius of the star exceeds its 
corotational radius (which is equivalent to the condition < 
PcdUsd), see Introduction) the linear velocity of the magnetosphere 
(which is assumed to co-rotate with the star) at the radius 
exceeds the sound speed in the surrounding material (the cor- 
responding state is usu ally referred to as supersonic propeller, 
lOavies & Pringlelll98lh . The heating rate of the atmosphere due 
to the propeller action by the star in this case significantly exceeds 
the rate of plasma cooling due to the bremsstrahlung emission and 
turbulent motions (for a discussion see also. Lamb et al. 1977, and 
references therein). The atmosphere during this state remains hot, 
T{r) Tf{{r). The rotational energy loss by the neutron star is con- 
vected up through the atmosphere by the turbulent motions and lost 
through its outer boundary. The spin-down rate of the star in this 
sta te (see Table in .Davies et al . 1979) is the same as that evalua ted 
by|m arionov & SunvaevI ( 1 19751) and used bv lUrpin et al.l ( Il998h in 
their calculations. This indicates that the evolution of the spin and 
the magnetic field of a spherically accreting neutron star up to a 
moment when its p eriod reaches /'cd(4 d) can be treated within the 
results presented bv lUrDinetaT] l ll998l) . In particular, this suggests 
that the surface field of the star to the end of the supersonic pro- 
peller state is only by a factor of 3 weaker than the initial field. 

As the spin period exceeds Pcdiki), the linear velocity of the 
magnetosphere at the radius becomes smaller than the sound 
speed corresponding to the temperature Tff{r^^). However, the ef- 
fective acceleration applied to the hot gas moving along the curved 
field lines. 



GM 



V2(r„) 



(9) 



is directed outwards from the star as long as the temperature of 
the material situated over the boundary is ( lElsner&Lamblll976l : 
p*Lrons & Lea 1976) 



rp(r„)>r„-o.irff(r„,). 



(10) 



Here r^„„ is the curvature radius of the field lines, 9 is the angle 
between the radius vector and the outward normal to the magneto- 
spheric boundary and VTi(''m) is the ion thermal velocity in the ma- 
terial located at the inner edge of the atmosphere. Furthermore, the 
magnetospheric boundary under the condition ilOt is stable with 
respect to interchange instabilities and the plasma penetration rate 
into the magnetospher e is limited to the rate of Bohm diffusion 
l Elsner& Lambl'l984'). which under the conditions of interest is 
(Ikhsa nov & Pustilnik 1996) 



TOb < 10" gs"' Q-, 



-1 1/2 1/7 



0.3S„ 



-1/14 



TO 



1015 



gs 



11/14 



(11) 



Here tto i = tt/0.1 is the diffus ion efficiency, which is normalized 
following iGosling etal] l ll99lh . Thus, the condition r,„ < ap- 
pears to be necessary, but not sufficient for a direct accretion with 
the rate of ~ TO onto the stellar surface to start. In addition, it is 
required that the cooling of plasma at the base of the atmosphere is 

more effective tha n the heating. 

As shown bv lOavies & Pringld ( Il98lh . the cooling of the at- 
mospheric plasma is governed by the bremsstrahlung radiation and 
the convective motion. For these processes to dominate the energy 
input into the atmosphere due to the propeller action by the star, the 
spin period of the star s hould be P^ <: Ph i, where P^r is the so-called 
break period, which is ( llkhsanovll200T3) 



PUki)^ 2000 m 



-4/21 



g(fsd) 

0.3S„ 



16/21 



TO 



1015 gs- 



-5/7 



(12) 



Under the conditions of interest, the break period significantly ex- 
ceeds f cd('sd)- This means that a spherically accreting neutron star 
is able to switch its state from the supersonic propeller to accretor 
only via an additional intermediate state, which is called the sub- 
sonic propeller. The star during this state remains surrounded by 
the hot adiabatic atmosp here and its spin period i ncreases from P^^ 
to Pbr on a time scale of JPavies &Pringlelll98ll) 



Td ^ 4 X 10^^ Us m-'^ 



0.3Ber 



-8/7 



TO 



'gs 



-3/7 



yr. 



(13) 



The accretion luminosity of the star during this state is limited to < 
TOgGM/rns, which under the conditions of interest is < 10"" erg s"i . 
Finally, the surface magnetic field of the star in the subsonic pro- 
peller state would not evolve significantly since Td « fsd- 



4 LONG-PERIOD PULSARS 

The results of incorporation of th e subsonic p r opelle r state into the 
evolutionary tracks presented bv lUrpin et al.l i l 19981) are shown in 
Tab[2] (the corresponding evolutionary scenario will be further re- 
ferred to as DFP-scenario). The second column of the Table shows 
the ratio of the average (on a time scale of f > Td) strength of the 
stellar wind during the spin-down epoch, < TO >, to the average 
mass capture rate by the neutron star inferred from the X-ray lumi- 
nosity, < TOx >= r„s < Lx > /GM. The third column represents 
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the ratio of the initial surface field of the neutron star to the critical 
field. 

The long-period pulsars in Table[2] are divided into three 
groups. The first group contains 4 sources whose origin within 
DFP-scenario can be explained provided the average strength of 
the wind during the system evolution has not been changed signif- 
icantly and that the initial magnetic field of the neutron stars was 
subcritical. The value of the ratio Bg/B^y for these sources has been 
evaluated from Eq. ( 124b by setting < SO? > / < 50?^ >= 1 and taking 
into account that fi(fsd) - 0.3So- 

An application of DFP-scenario to J0103. 6-7201 and 
J0146.9+6121 (the second group) leads us to a conclusion that ei- 
ther the wind during a previous epoch was by a factor of 2-5 weaker 
than that inferred from the X-ray luminosity of these sources or the 
value of 6(fsd) wa s a factor of 2 larger than the average one es- 
timated bv lUrpin et al. (1998), or both. According to th e results 
reported bv I Vink et aL t200a) and iFullerton et al.1 BoO^ it is not 
unusual for massive stars to lose material with a higher rate at later 
stages of their evolution. It also cannot be excluded that the X-ray 
luminosity of these sources is slightly overestimated or has been 
measured during a high state of activity of the massive compan- 
ions, which usually lasts 3-20 years. In this light, the assumption 
about a weaker wind in a previous epoch of 10103.6-7201 and 
J0146.9+6121 seems to be rather reasonable. On the other hand, 
the duration of the spin-down epoch of neutron stars with the ini- 
tial field ~ Ba in the wind of ~ 10'^g s~' is by a factor of 2-3 
shorter than the average one evaluated bv lUrpin et alj ( [l998l) . This 
indicates that the assumption about a higher value of fi(4d) has cur- 
tain theoretical grounds as well. However, both of these arguments 
are not compelling to exclude a possibility that the neutron stars in 
J0103. 6-7201 and J0146.9+6121 were born as magnetars. 

The third 'group' is represented by only one (but the most puz- 
zling) source 2S 01 14-1-650. The origin of this source within DFP- 
scenario can be explained provided the strength of the wind in a 
previous epoch was by a factor of ~ 30 smaller than the mass cap- 
ture rate by the neutron star inferred from the observed X-ray lumi- 
nosity. Indeed, for the condition Pbr = lO"* s in the case Bq < Scr to 
be satisfied the strength of the wind should be 3j? < 10'* g s"' . The 
neutron star under these conditions is able to reach a period of 10"* s 
within a time fsd < fms (see Fig. 2 in Urpin et al. 1998, and Ea.ll3t. 

On the ot her hand, an attempt to interp ret this source within 
the scenario of iLi & van den Heuvell il99% leads us to a signifi- 
cantly stronger limitation to the value of SJi. Indeed, the neutron 
star within their scenario is presumed to switch to the accretor state 
as soon as its spin period reaches /'cd(rsd)- Setting Bq ~ 10'^ G 
they have evaluated the spin-down time scale as tsd < 10' yr 
and, correspondingly, B(tsd) ~ (1 - 4) x lO'^'G, which implies 
/^(^sd) ~ 10^^^ yU32Gcm''. Putting this value to Eq. (TJ and solv- 
ing the equation for W((t^^) one finds 9jt(?sd) S 2 x 10''' gs"' x 

K'^^rrr^^^n\2 (^'cd(fsd)/ lO"* s) . This means that the interpretation 
of the long spin period of the neutron star in 2S 01 14+650 within 
the magnetar hypothesis requires the strength of the wind in a pre- 
vious epoch to be by a factor of 250 smaller than that inferred from 
the X-ray luminosity of the source. 

One of the ways to improve the situation is to assume that 
the evolutionary tracks of magnetars contain the subsonic pro- 
peller state. In this case the neutron star can comfortably reach 
a period of 10"* s even if the strength of the wind is ~ 10" gs"' 
(see Sect.O. However, a necessity to invoke this assumption is 
not obvious. Accor ding to the evolutionary tracks reported by 
iMevnet et all ( Il994h . a possible progenitor of the Bl-supergiant 



Table 2. Parameters of long-period pulsars inferred from DFP-scenario 



Name 


< jji > / < iUix > 




J170006-4157 


1 


0.08 


0352+309/X Per 


1 


0.3 


JI037.5-5647 


1 


0.3 


J2239.3+6116 


1 


0.1 


J0103.6-7201 


0.2-1 


1 


J0146.9+6121/V831Cas 


0.6 


1 


2S 0114+650* 


0.03 


1 



* See, however, discussion in the text. 

was a 09.5 V star. The typic al mass-loss rat e of this type stars (see, 
e.g.. lCrowther. Lennon & Wa lborn 2006; Fullerton et al.ll2006l. and 
references therein) is by a factor of 10-50 smaller than the mass- 
loss rate of the massive companion in 2S 01 14+650 inferred from 
the optical observations dAab & Bvchkovalll984ljReig et al.lll996h . 
This indicates that the value of the ratio < 9)1 > / < SJlx >, given in 
Tab.[2]has good observational grounds and, therefore, allows us to 
interpret this object in terms of DFP-scenario assuming the initial 
magnetic field of the neutron star to be subcritical. Finally, I would 
like to remind that the association of the 10"* s pulsations with the 
spin period of the neutron star cannot be considered as finally jus- 
tified (see Sect.|2j. If the spin period of the neutron star is indeed 
smaller than that currently adopted the origin of this source can be 
interpreted within DFP-scenario in terms of the subcritical initial 
field almost without any additional assumptions. 



5 DISCUSSION 

The main result of this paper is that the problem posed by the ex- 
istence of relatively bright persistent long-period X-ray pulsars can 
be solved by an incorporation of the subsonic propeller state into 
the evolutio nary tracks of neut ron stars in massive close binaries 
computed bv lUrpin et alj ( fl99^ . This solution implies that the neu- 
tron stars in these systems undergo spherical accretion and their 
initial magnetic field is subcritical. The spin periods of the neu- 
tron stars in this case are close to the equilibrium period expressed 
by Eq. ([8j provided the relative velocity between the star and sur- 
rounding gas is Vrei ~ 400 - 800 km s"'. The only exception is the 
neutron star in 2S 01 14+650 whose equilibrium period for a rea- 
sonable estimate of V,ei is by an order of magnitude smaller than 
10* s. It should be, however, noted that the association of the 10* s 
pulsations with the spin period of the neutron star in this system 
remains so far controversial. 

The evolutionary scenario presented seriously weakens the hy- 
pothesis that the neutron stars in the long period pulsars are magne- 
tar descendants. In particular, it is shown that the interpretation of 
the long-period pulsar 2S 01 14+650 within this hypothesis is not 
possible unless an assumption about a significantly (by a factor 
of 250) weaker wind in a previous epoch is invoked. Furthermore, 
this hypothesis implies that magnetars form during the first super- 
nova explosion in a massive binary system, which is not consistent 
with the current views on the process of magnetar formation. At the 
same time the above arguments are not compelling to exclude the 
magnetar hypothesis completely and therefore, a brief discussion 
about possible appearances of magnetars in close binaries appears 
to be rather reasonable. 
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The main evolutionary stages of a magnetar witli an initial 
magnetic field of fig ~ 2 x lO'"" G (jUo = 10'^ Gcm^) in a binary 
system are as follows (for a discussion see, e.g., IlkhsanovlliooTbL 
and references therein). The ejector stage, which lasts 



Ta~ IOV33 '45 3jf75" V, yr, 

and ends as the spin period of the neutron star is 



53 wql^' Vg s 



-1/4 , 



(14) 



(15) 



where Vg = V,-ei/10'* cms" . 

The parameters of the propeller stage depend on the geome- 
try of the accretion flow beyond the neutron star's magnetospheric 
boundary. If the magnetospheric radius of the star is smaller than 
the so called circularization radius, which is defined as 



J' 



where 



(16) 



(17) 



is the rate of accretion of angular momentum, the formation of a 
disc would be expected. Combining Eqs. J16b and ( I17t and solving 
the result for V^ei one finds the condition for a disc formation as 



v., < v.- - 75 ^(f ^-32'" ^^5o' »J!f kms-'. 



(18) 



It is taken into account here that the magnetic field of a magnetar on 
a time sca le of ~ 10^ yr dec reases by almost an order of magnitude 
(see, e.g.. lColDi et al.ll200Cli) . If this condition is satisfied the spin 
period of the neutron star during the propeller stage decreases to a 
value of 



on a time scale of 



■3/2 ,6/7 ,„-5/7 s_,jj-3/7 



(19) 



\l/2 

-1 ) y-' 



(20) 



r-. 2x10^ (^) 

and the magnetar switches its state to an accretor. Here K0.5 = k/0.5. 

If the plasma flow beyond the magnetospheric boundary has 
a spherical geometry the duration of the supersonic propeller state 
would be 

T<^P> - 6 X 10^ I45 ti-^i V-^"^ yr. (21) 

The spin period of the star during this time increases to a value of 

< p(f ^ 103 m-'" m;l" s. (22) 

As the star switches to the subsonic propeller state its spin period 
increases on a time scale of 



Td ~ 200 yUjj''' hs m^''' ffll^j" yr. 



to a value of 



A ,6/7 ^-5/7 „ -4/21 



1.1 X IOV32 ™ 



IS 



(23) 



(24) 



The above estimates suggest that the spin period of magnetars 
accreting material from a disc unlikely exceeds 500Ljj'' s, where 
L35 is the X-ray luminosity of the source expressed in units of 
10^^ erg s"' . If the time scale of the magnetic field decay of magne- 
tars exceeds 3 x lO"* yr a probability of observational identification 
of these stars in the magnetar stage is not zero. On the other hand, 
it appears to be rather small as the magnetic field of an accreting 
star decays m ore rapidly than th at of an isolated one (for a discus- 
sion, see, e.g.. lUrDin et alj[l998l and references therein). The spin 



period of magnetars in the accretor state evolves to the equilibrium 
period expressed by Eq. (O and the time scale of its evolution is 
determined by the decay rate of the stellar magnetic field. 

The magnetars undergoing spherical accretion can be spun- 
down to a period of ~ 2 x 10''L35'' s. However, the spin-down time 
scale in this case proves to be comparable with the upper limits 
to the characteristic time of the magnetic field decay. Therefore, a 
possibility to observe an accretion-powered magnetar undergoing 
spherical accretion is almost negligible. 

Finally, the appearance of close binaries containing magnetar 
descendants of an age ~ /„,s would be similar to that of close bi- 
naries in which the initial field of the neutron star was under the 
critical value. In contrast, the evolutionary tracks and appearance 
of neutron stars undergoing disc and spherical accretion differ sig- 
nificantly. 
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